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Abstract. We present some results obtained in collaboration with prof. Piero D'Ancona 
concerning global existence for the 3D cubic non linear massless Dirac equation with a 
potential for small initial data in H 1 with slight additional assumptions. The new crucial 
tool is given by the proof of some refined endpoint Strichartz estimates. 



1. Introduction 

We consider the 3D massless Dirac equation perturbed with a small potential V 



(1.1) 



\iu t + Vu + Vu = P 3 (u) 
\u(0) = f{x) £ H s , 

where the Dirac operator is defined as 
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V = l/i^2 ®kdk = -i(a • V) 
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the nonlinear term Ps(u) = \u\ 3 and V is a 4 x 4 hermitian matrix. Notice that 
(1.2) aja k + a k aj = 25 jk l, j,k = 1,2,3; 

The unperturbed nonlinear Dirac equation is important in relativistic quantum mechanics, 
and was studied in a number of works (see e.g. [13], [3], [12], [11], [4], [10], [9]). In 
particular, it is well known that the cubic nonlinearity is critical 1 for solvability in the 
energy space H 1 ; global existence in H 1 is still an open problem even for small initial 
data, while the case of subcritical spaces H s , s > 1 was settled in the positive in [4], [10]. 
The standard tool when proving global wellposedness for general dispersive equations in a 



^The general framework when dealing with low regularity well posedness for nonlinear dispersive equa- 
tions is to be thought as follows: the scaling of the operator, the functional setting and the structure of the 
nonlinearity yield a threshold regularity such that if the initial data has more regularity than this (i.e. is 
subcritical) it is possible to prove global wellposedness with standard contraction arguments. If the initial 
data is instead supercritical then we have various degrees of illposedness, while the critical case is typically 
harder and many different things can happen. 

1 



2 



FEDERICO CACCIAFESTA 



low regularity setting, at least in the subcritical case, relies on some space-time estimates 
(Strichartz estimates) that allow the application of the classical fixed point method on 
some suitable functional space; the critical case may be corresponding to the endpoint 
estimate and thus may require a different approach. 

To begin with, let us give a look to the free homogeneous case. Notice that property 

(1.2) yields V 2 = -J 4 A, so that 

(1.3) (id t -V)(id t + V)u = -Uu. 

This means that every solution to the free Dirac equation is a solution to a vectorial wave 
equation with suitable initial conditions. Despite the rich structure of T> it is thus easy, 
relying on the well known theory for the free wave equation, to write the Dirac linear 
propagator as 

(1-4) ^f = cos(t\D\)f + i^^-Vf. 



and the Strichartz estimates 



that hold for every couple (p, q) such that 

2 2 

- + - = 1, 2<p<oo, oo>g>2. 
P q 

We recall that estimates (1.5) capture two different informations: 

Locally in time they describe a type of smoothing effect, reflected in a gain of integra- 
bility with respect to some suitable norm. 

Globally in time they describe a decay effect, meaning that some spacial norm of a so- 
lution must decay to zero as t — > oo, at least in some LP time-averaged sense. 

The choice of the couple (p, q) = (2, oo) yields the so called endpoint estimate 
(1-6) \\e lW f\\ L ^ < II/Ha! 

that is known to fail, as the corresponding one for the 3D wave flow (see [8]). If (1.6) 
were true one could prove global wellposedness for problem (1.1) (at least for V = 0) with 
standard methods. Considering the solution map 

v i— > = e iW f + i f e^-^Psivit'^dt' 

one would have indeed 



o 



-t 

e 



10 



L 2 L° 



/+oo 
\\e lW e- u ' v P 3 (v(t'))\\ L2L ^dt'<\\v 3 \\ L , m 
-oo 



THE CUBIC NONLINEAR DIRAC EQUATION 



3 



that in conjuction with the conservation of H norm would imply 

\Mv)\\ LrH i + ||*(u)|| £ a L « < WfWm + l|«IU«Hi||«lli2£cc. 

In other words, a contraction argument in the norm || • || £ 2 £ oo + || • || £ oo^i would be enough 
to prove global existence of small H l solutions to (1.1). 

The failure of estimate (1.6) prevents us from the application of this machinery. Our aim 
here will thus be to provide some "weak endpoint" estimates, i.e. some suitable refinements 
of estimate (1.6) in view of obtaining global wellposedness for problem (1.1) in the critical 
case with some slight additional hypothesis on the initial data. 

We shall here present two approaches to the problem, that in both cases will prove, at 
least, global existence for problem (1.1) for radial H 1 initial data. 



2. An "algebraic" approach 

In [8] the authors noticed that the endpoint estimate for the 3D wave flow can be 
recovered if the initial data are assumed to be radial. As it can be easily seen indeed, for 
a radial / we have, for all x, 

\ D \ M J\\x\-t\ 

where M(g) is the Hardy maximal function of the function g(s) = sf(s), so that 

<M(g)(t) 

and thus by a standard maximal estimate 

^ IIS'IIl^R = \\sf(s)\\ L 2 {R ^ ||/||X2(R3) 

(analogous arguments hold in dimension n > 3, see [5]). This remark suggests the chance 
of slightly improving the range of admissible exponents when requiring some additional 
structure on the initial data. Even if the Dirac operator does not preserve the radiality of 
functions, we are indeed able to prove the following 

Proposition 2.1. Let n = 3 ad let f belong to the space 

(2.1) ii 1 = {/i + Vf 2 , h e i^OR 3 ), h G i? 2 (M 3 ), fi, h radial}. 

Then the following endpoint Strichartz estimate holds: 

(2-2) V W f\\LtL~ < WfWm- 

Our proof of this result, already contained in [9] and mainly relying on Fourier transform 
in radial coordinates, takes the remarkable advantage of being adaptable also to handle the 
non homogenous term. With minor modifications we are in fact able to prove the following 
mixed Strichartz-smoothing estimate 



sin(i|L>|) 
\D\ 



sin(t|D|) 
\D\ 
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Proposition 2.2. Let n = 3 and assume F(t,x) has the structure 
(2.3) F(t,x) = F 1 (\x\)I 4 + i(a-x)F 2 (\x\). 

Then the following estimate holds 



(2.4) 



f e i{t - s)v F{s) ds 
Jo 



< 



(x)2~ 



\D\F\ 



This result is crucial in view of obtaining endpoint Strichartz estimates for the potential- 
perturbed Dirac flow, and to keep the loss of derivatives less or equal than 1. Propositions 
2.1 and 2.2 yield in fact the following Theorem, a proof of which will be also sketched. 

Theorem 2.3 (Endpoint estimate for the perturbed Dirac flow). Let V{x) be a 4 x 4 
matrix of the form 

(2.5) V(x) = V 1 (\x\)I 4 + ip(a-x)V 2 (\x\), Vi,y 2 :R + ^M. 

Assume that for some a > 1 and some sufficiently small 5 > 

5 6 



(2.6) 



\V(x)\ < 



(x)y 2 +wl' 2 ' 

with w a {x) 1 / 2 = \x\{\ + |log|x||) ff . 

Then the following endpoint Strichartz estimate 



\W(x)\ < 



(2.7) 



D it(V+V) 



/lb 



H 1 



holds for all initial data f £ rl 1 fl H . 



Proof. By Duhamel's formula we write the integral form 

r-t 

MV+V) f = e itV f + i I J(ts)T> 



u 



(Vu)ds 







so that we have 
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Jt(V+V) 



3 i(t-s)X> 



(V(s)e 



is{V+V) 



f)ds 



Applying Propositions 2.1 and 2.2 we can continue estimating with 

<ll/bi + ll^)^I^K^)llxfx- 
From the L 2 boundedness of the Riesz operator |Z)| -1 V we have 



< 



-(yv)u\ 



+ \\{x)2+V(Vu) 



\L\Ll- 



To conclude we need the following smoothing estimates (the first of which is proved in 
D'Ancona-Fanelli '08, the second directly follows), which hold with more general assump- 
tions on the potential V: 



L 2 ; 
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. ||^ 1/2 Ve^)/|| L?L , < \\f\\ Hl . 
Thus multiplying and dividing by 

w a {x) 1 ' 2 = |x|(l + |log|x||) CT 

in the smoothing terms yields 

. \\{x)^{VV)u\\ L , Ll < \\{x)^+ W l /2 VV\\ L ~ ■ \\ W ? /2 u\\ L * Ll 

<\\fh^ 

. \\(x)^V(Vu)\\ LUl < \\(x)V 2 +wl /2 V\\ L ~ ■ \\w^ /2 Vu\\ L i Ll 

< II/IIai 

and this concludes the proof. □ 

Remark 2.1. Notice that hypothesis (2.6), that is needed in view of applying proposition 
2.2, is fairly natural: the structure required describes in fact electric potentials and particles 
with anomalous magnetic momentum. 

It would be now tempting to apply this result to prove global existence for problem (1.1) 
for small initial data in ^PiH 1 with standard fixed-point techniques; since such a method 
is iterative, what we would need is that class to be invariant under the action of the cubic 
non linearity Ps(u), but this fact is unfortunately not true even with the simplest choice 
P 3 (u) = ((3u,u)u. 

Our next step will thus be to show the existence of some proper subspace of the set 
n H 1 on which the following conditions are satisfied: 

• The operator T> + V is well-defined and selfadjoint, 

• The action of Ps(u) is invariant. 

The classical theory of partial wave subspaces (see [15]) allows in fact, analyzing the 
structure of the Dirac operator in radial coordinates and using spherical harmonics, to 
decompose the space L 2 (M 3 ) 4 in a direct sum of Hilbert spaces that are left invariant 
under the action of the potential perturbed Dirac operator. Moreover, an explicit insight 
of these spaces shows that "some of them" are left invariant by the action of the standard 
cubic non linearity as well. We collect all this results in the following proposition, referring 
to [15] and [1] for notation, details of the statement and a proof. 

Proposition 2.4. There exists a family of 2- dimensional Hilbert spaces T-L mj ,kj (the "par- 
tial wave subspaces") such that 

L 2 (R 3 ) 4 0L 2 ((O,oo),dr) ®U mj , kj 

and the operator D + V , with V of the form (2.5) acts and is self-adjoint on such spaces. 

Moreover for j = 1/2 the spaces Hm 1 i 2 ,k 1 /2 are ^ e ft invariant a ^ so by the action of the 
cubic nonlinearity P-s(u) = (f3u,u)u. 
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Remark 2.2. Each of these partial wave subspaces has a basis {^^. ki^ m . that can 
be explicitly written using Legendre polynomials, satisfying the property that 



The action of the operator T> + V can thus be easily represented in coordinates with respect 
to such basis, and in this contest is commonly referred to as the radial Dirac operator. 

Remark 2.3. The second part of the statement is justified by the fact that for j = 1/2 it 
is easy to evaluate the action of P3 on the functions k . , <3?~ k . , that in this particular 
case have a fairly simple structure. To give a clearer picture of the situation, consider a.e. 
the triple 2 (j, rrij, kj) = (1/2, 1/2, 1); then we have 



/ 



1/2,1 



2v^F 
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As it is easily seen, this yields the fact that for a generic function u € L 2 ((0, 00), dr)®%\/2^ 
that will thus be written as 

u(r, 9, cf>) = u + (r)$+ /2)1 (0, <f>) + u~(r)^ /21 (9, </>) 
for some radial functions u + , u~, the nonlinear term reads as 



(/3u,u) 



1 

47T 
1 



cos 2 8 u + (r) 2 — — sin 2 6 u + (r) 2 — —u (r) 2 
4tt Air 



-(u+{rf+u-{rf) 



and so has no angular component. This proves that the action of the standard cubic 
nonlinearity is invariant on the space %i/2,i- 

This last remark provides the final tool we needed to state our global existence result, 
that now can be proved with completely standard techniques. 

Theorem 2.5. Let P%(u) = {/3u,u)u and the potential V satisfying (2.5)- (2.6). Then for 
every initial data f € i7 1 ((0, 00), dr) ® Hm 1 / 2 ,ki/2> w ith sufficiently small H 1 norm, there 
exists a unique global solution u(t,x) to problem (1.1) in the class Ct(R, H 1 ) n L 2 (R, L°°) . 



Notice that for a fixed j = jr, 5, ... there will be 4j + 2 possible choices of rrij, kj. More precisely the 
range is rrij = -j, + +j and k 3 = -(j + 1/2), +(j + 1/2). 



THE CUBIC NONLINEAR DIRAC EQUATION 



7 



3. Angular regularity estimates 

In recent years many papers (see e.g. [14], [6], [7], [9] and references therein) have been 
devoted to the improvement of Strichartz estimates using angular regularity. In particular 
in [9] it is proved the following estimate for the wave propagator 

(3-1) n = 3, \\j m f\\ L 2 L? > L p < y/p-\\\D\f\\ L 2, Vp<oo 

where we are using the natural notation 



and 



L?L b u =sup||/(r • )|| £ 6(s»-i). 

r>0 



Notice that the norm at the left hand side distinguishes between the integrability in the 
radial and tangential directions. Using estimate (3.1), Machihara et al. were able to prove 
global well posedness for problem (1.1) with V = for small .fp-norm data with slight 
additional angular regularity, and in particular for all radial H 1 data. This is especially 
interesting since, as pointed out in the introduction, radial data do not correspond to radial 
solution for the Dirac equation (due to the fact that the operator D does not commute 
with rotations of K 3 ). 

Estimate (3.1) gives a bound for the standard L 2 L°° norm via Sobolev embedding on 
the unit sphere § 2 

(3.2) V m f\\ L *L~ < WK^fW^US, £ WimUh*, P>\ 

where the angular derivative operator is defined in terms of the Laplace-Beltrami op- 
erator on W 1 ' 1 as A* = (1 - A§„-i) s / 2 . 

Our main goal here is to extend this group of results to the Dirac equation perturbed 
with a small potential V(x). We consider first the linear equation 

(3.3) iut = Vu + V(x)u + F(t,x). 

The perturbative term Vu can not be handled using the inhomogeneous version of (3.1) 
because of the loss of derivatives. Instead, we shall need to prove new mixed Strichartz- 
smoothing estimates. 

We collect in the following Theorem the Strichartz estimates we are able to prove both 
for the wave and Dirac equations (see [2]). 



Theorem 3.1. Let n > 3 and 

(3.4) a, 



10 if n = 3 
\l-f i/n>4. 

Then for every s > the following estimates for the free wave equation hold 
(3-5) \\K<& m f\\a!L-i* , < l|A* +CT "/|l 



H~ 
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(3.6) 



Af 



M-s)\d\ f ^ x ^ 



< 



(X 2 



^■D\^A S +^F\ 



J 2 r cxi r 2 



77ie corresponding estimates for the 3D free Dirac equation are 
(3-7) \\Ke UV f\\ L ^ Li < WKfWfr, 



(3.8) 



A , f l e i(t-t')v F ( t ' :X ) dt < 
Jo 



< \\(x)-i + \D\KF\\ L i L2x . 



Remark 3.1. As a byproduct of our proof we have thus obtained the following endpoint 
estimates for the wave flow with gain of angular regularity: 

(3-9) n > 3, \\e UlDl f\\ L 2 L?aLl < 11^-/11^ 

where a n is as in (3.4). Although this was not the main purpose here, it is interesting to 
compare (3.9) with known results. In dimension n = 3, estimate (3.9) is just a special case 
of Theorem l.l-III in [9] where (3.9) is proved with 03 = — |; it is not known if this value 
is sharp, however in the same paper it is proved that the estimate is false for 03 < — |. On 
the other hand, to our knowledge, estimate (3.9) for n > 4 and (3.6) for n > 3 are new. 
The literature on these kind of estimates is extensive and we refer to [6] , and the references 
therein for further information. 

The next step is to prove suitable smoothing estimates for the Dirac equation with 
potential. By a perturbative argument we obtain the following endpoint estimates for the 
linear flows: 

Theorem 3.2. Assume that the hermitian matrix V{x) satisfies, for 5 sufficiently small, 
C arbitrary and a > 1, with v(x) = \x\ 2 1 log \x\\ 2 + + 

(3.10) |v(,)| < J^, |w(,)| < ^r y 

Then the perturbed Dirac flow satisfies the endpoint Strichartz estimate 

(3.11) \w tiv+v) f\y t L?Li< ii/iifli. 

If the potential satisfies the stronger assumptions: for some s > 1, 

(3.12) ||A* • )|| L 2 (§ 2) < -J-,, ||A« W(|x| • )|| L 2 (§ 2) < ^- 
then we have the endpoint estimate with angular regularity 

(3-13) \\Ke lt[V+V) f\\ L ^Ll < WKfWm 

and the energy estimate with angular regularity 
(3-14) \\A s u e it ^f\\ Lrm < \\KfWm 
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Remark 3.2. Notice the difference in the assumptions on V between Theorems 2.3 and 3.2: 
the hypothesis on structure (2.5) here is substituted by (3.12) which ensures some minimal 
angular regularity needed. 

We can finally apply Theorem 3.2 to the nonlinear equation (1.1) to obtain: 

Theorem 3.3. Consider the perturbed Dirac system (1.1), where the 4x4 matrix valued 
potential V{x) is hermitian and satisfies assumptions (3.12). Let P^(u,u) be a C 4 -valued 
homogeneous cubic polynomial. Then for any s > 1 there exists eo such that for all initial 
data satisfying 

(3-15) WKfWm < eo 

the Cauchy problem (1.1) admits a unique global solution u £ CH 1 n L 2 L°° with K s u u £ 
L^H 1 . 

In particular, problem (1.1) has a global unique solution for all radial data with suffi- 
ciently small H 1 norm. 

Remark 3.3. It is clear that our methods can also be applied to nonlinear wave equations 
perturbed with potentials, and allow to prove global well posedness for some types of 
critical nonlinearities. 

Remark 3.4. We did not strive for the sharpest condition on the potential V, which can be 
improved at the price of additional technicalities which we prefer to skip here. Moreover, 
differently from the previous chapter in which the structure of the non linear term was 
essential, the result can be extended to more general cubic nonlinearities |P3(u)| ~ | ix| 3 . 

Remark 3.5. Notice that we need an angular regularity s > 1 on the data, higher than the 
s > assumed in the result of [9]. It is possible to relax our assumptions to s > 0; the 
only additional tool we would need to prove is a Moser-type product estimate 

I|A*M|| L 2 < lMMI A >lk + l|A>|| L 2 ||t;|| L oo, s > o 

and an analogous one for A^DKuv). This would require a fair amount of calculus on the 
sphere S 2 , and here we preferred to use the conceptually much simpler algebra property of 
H 3 ^- 1 ) for s > 2fi. 
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